Electrons and Phonons in YbCg 
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The electronic structure and selected zone-center phonons in Yb graphite intercalation compound 
(YbCe) are investigated theoretically using density functional calculations and LDA+U approxima- 
tion for Coulomb correlations in the / shell, and experimentally using angle-resolved photoemission. 
We find that both in LDA and LDA+U approach the Yb / states are fully occupied providing no 
evidence for mixed-valent behavior. The obtained theoretical results are in good agreement with 
photoemission experiments. The 4/ states are considerably hybridized both with the Yb 5d and 
C 2p states resulting in a substantial admixture of Yb / at the Fermi level. Soft Yb phonons, 
given a noticeable presence of the Yb states at the Fermi level, are probably responsible for the 
superconductivity recently reported in YbCe. 



For a long time graphite intercalation compounds 
(GICs) have been attracting substantial attention due to 
their layered quasi-two dimensional (2D) structure and 
the resulting large anisotropy of their electric and elec- 
tronic properties. 1 ' 2 Some of GICs are known to be su- 
perconducting at low temperatures. 3 A renewed interest 
in superconducting phenomena in intercalated graphites 
is due to the recent discoveries of superconductivity in 
MgB2 and B X C, electronically related to graphite. 4 An 
interesting example of graphite intercalation compounds 
is YbCg. First, it is one of the only four known bulk 
intercalations with 4/ metals: YbC6, EuC6, SmCg, and 
TmCg which may indicate that / electrons play certain 
role in bonding in these compounds. Likewise, YCg has 
not been synthesized so far, although Y behaves simi- 
lar to lanthanides in systems with localized / electrons 
(cf. high-T c cuprates). Last but not least, very recently 
superconductivity was discovered in YbC6- 6 

Yb is known to form mixed-valent compounds, 7 so one 
may think mixed-valent physics is operative in YbCg as 
well. Therefore, it is of interest to investigate the elec- 
tronic structure of YbCg numerically and experimentally. 
Should a good agreement between the theory and ex- 
periment be established, the former can be used to gain 
some insight into possible mechanism for superconductiv- 
ity. Experimentally, Yb valence in the above compound 
can be obtained by recording the Yb 4/ photoemission 
(PE) spectra. Due to the large Coulomb-correlation en- 
ergy, the energy positions of the PE signals for 4/ con- 
figurations with different electron occupations are shifted 
by several eV with respect to each other. Therefore, the 
contributions to the PE from different 4/ configurations 
can easily be discriminated. 8 From the intensity ratio of 
these contributions the information about the valence can 
be derived. 

In this paper we present ab initio calculations of the 



electronic structure of YbC6 using both local density ap- 
proximation (LDA) and LDA+U approach, which ac- 
counts for Coulomb correlations inside the / shell. The 
calculated energies for the Yb 4/ states are in good agree- 
ment with the results of photoemission experiments giv- 
ing no evidence for the mixed-valent behavior. We also 
report first principles calculations of selected zone center 
phonons. We find soft Yb-derived phonons to be likely 
responsible for the superconductivity reported in YbCg. 

For the calculations, we used experimental crystal 
structure P6/mmc, with C occupying 12i and Yb 2d 
positions, with the lattice parameters a — 4.32 A and 
c = 9.1 A (Fig. 1). A full potential linear augmented 
plane wave method (LAPW) was used 9 with the follow- 
ing setup: The APW sphere radii were taken as 2.5 and 
1.3 bohr, the cutoff parameter i?A" max = 7, and local 
orbitals were used for Yb s and p and for C s, to re- 
duce the linearization error, and to improve convergence 
in RK max . Linear tetrahedron method was employed for 
the Brillouin zone (BZ) integration, with the k mesh up 
to 11x11 x4 divisions. Spin-orbit coupling was included 
on the second variational basis. 9 The Ceperley- Alder 
exchange-correlation potential was used in the LDA part 
of the calculations. Finally, Hubbard correlations in the 
/ shell were taken into account using the fully localized 10 
version of the LDA+U formalism, with the parameters 
U = 0.4 Ry and J = 0.05 Ry, estimated by modifying 
occupation numbers in the quasiatomic loop in a LMTO 
(linear muffin-tin orbital) program. 

The Yb-GIC samples for the PE experiments were 
prepared in situ by thermally driven surface reaction of 
deposited Yb overlayers with graphite (0001) substrate 
as described elsewhere. 11 The measurements were per- 
formed using synchrotron radiation from various beam- 
lines at BESSY/Bcrlin electron storage ring. Angle- 
resolved photoemission spectra (ARPES) were taken 
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with 33 eV photon energy using rotative hemispherical 
electron energy analyzers (VSW- ARIES) tuned to an en- 
ergy resolution 150 meV (full width at half maximum, 
FWHM) and an angular resolution of 1°. Basic pressure 
during measurements was always better than lxl0~ 10 
mbar. 




FIG. 1. Crystal structure of YbCe- Large spheres denote 
Yb and small spheres denote C. (color online) 

The band structure calculated within the LDA+U 
scheme is shown in Fig. 2. Probably the most intriguing 
observation is that all / states appear to be fully occu- 
pied, giving evidence to firmly divalent Yb. This seems to 
be related to the fact that straight LDA calculation (not 
shown), without any account of intra- /-shell Hubbard- 
type correlations, still place all / bands below the Fermi 
level, although, of course, much higher in energy than in 
LDA+U. 
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FIG. 2. LDA+U band structure of YbC 6 . The left panel 
shows the partial Yb-d character, and the right panel Yb-/ 
character, (color online) 



The spin-orbit interaction splits the / states into two 
manifolds, located 1.10 and 2.30 eV below the Fermi 
level (Ep). The corresponding total electronic density of 
states (DOS) with two strong peaks originating in the 4/ 
states is demonstrated in Fig. 3. Calculated positions of 
the 4/ states in lanthanide systems depend very much on 
specific choice of U and J values, which are often used as 
fitting parameters to achieve better correspondence with 
the experimental data. It is very important to note that 
in the present study the U and J values were estimated 
from the first principles and no attempts were done to fit 
or tune this parameters. 
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FIG. 3. Electronic density of states of YbC6- The main 
graph shows the total DOS on a per-formula-unit basis and 
the insert shows, in the order from the lowest to highest DOS 
at E — Ef = 1 eV, the partial contribution of Yb /, Yb d, 
C p character inside the corresponding APW spheres and the 
contribution of interstitial states, (color online) 

Most importantly, the qualitative prediction of purely 
divalent Yb in Yb-GIC is unambiguously supported by 
the PE data. Fig. 4 depicts a series of angle-resolved 
PE spectra taken along the T-K'-M' direction in the 
BZ of graphite. All four characteristic occupied bands 
of graphite {ui v , (T2v, &3vi and 7Ti„) are clearly seen in 
the figure pointing to the high-quality intercalation com- 
pound under consideration. 11 The Fermi-energy peak ap- 
pearing in the region close to the if' point stems from 
the bottom of the 7Tq band of graphite, which becomes 
occupied in the GIC due to charge transfer of electrons 
from the intercalant atoms. The none-dispersive divalent 
contributions into the PE signal of the Yb 4/ electrons, 
which arc found at 1.14 eV and 2.41 eV binding energies 
(BEs) , are marked by two thin lines through the spectra. 
Would the divalent 4/ component of the final-state PE 
multiplet be observed at the Fermi energy, it will show 
that the final-state multiplet is energetically degenerate 
with the ground state that will provide a reason for ho- 
mogeneous mixed- valent behavior of the system. In the 
present case of the Yb-GIC, the divalent signal is ob- 
served too far from Ep to anticipate the mixed- valent 
properties of the compound. Correspondingly, no traces 
of a trivalent signal in the region of 5 to 12 eV BE, 8 which 
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one would expect for a mixed-valent system, are seen in 
the photoemission spectra. 

The experimentally obtained binding energies of the 
divalent 4/ contributions are in excellent agreement with 
the theoretical values. Not only energies, but also line- 
shapes of the DOS and PE spectra are in agreement 
with each other. In the right panel of Fig. 4 we com- 
pare the angle-resolved PE spectrum sampling electronic 
states in the region of the M point in the BZ of the Yb- 
GIC with the local DOS (spectral function) calculated 
around this point. The calculated local DOS was broad- 
ened with a Lorcntzian to account for finite lifetime ef- 
fects. Good correspondence between the theoretical and 
experimental data is evident in the figure. In fact, the 
agreement is better than one expects given the approx- 
imations used. Nevertheless, even though such level of 
accuracy is likely fortuitous, this clearly indicates that 
the underlying physics is correctly described by the cal- 
culations, and, importantly, that the derived value of U 
and J are in the right ballpark. 
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FIG. 4. (Left) angle-resolved PE spectra of the Yb-GIC 
taken along the V-K'-M' direction in the BZ of graphite 
(r-M-r-M in the BZ of the Yb-GIC) at different polar angles 
Q. (Right) Comparison between the local DOS in the vicinity 
of the M point and the angle-resolved PE spectrum sampling 
this region. 

We shall now discuss the calculated electronic struc- 
ture in more detail. Even in the vicinity of the Fermi 
level, the bands of YbCg are qualitatively different from 
that of YCg (Fig. 2), indicating that the hybridization of 
the / states with the itinerant electrons should play an 
important role in transport properties in this material. 
There are 12 bands crossing the Fermi level. However, 
spin-orbit splitting at E^ is very small, so that one can 
safely speak about 6 distinct Fermi surfaces, correspond- 
ing to bands 1, 3, 5, 7, 9, and 11, shown in Fig. 5. An 



easy way to make the correspondence between the Fermi 
surfaces and the bands plotted in Fig. 2 is to count them 
from K to T: the only band crossing the Fermi level be- 
tween K and M is band 1(2), while the state at approx- 
imately 0.2 eV at V is band 11(12). Bands 1, 5, 7, and 9 
are quasi 2D, while bands 3 and 11 have considerable dis- 
persion along the z direction. It is instructive to look at 
the element specific character of these bands. Band 11, 
forming elliptical pockets around T, has more Yb char- 
acter than C, and namely Yb d z 2_ r 2 character. Not sur- 
prisingly, this is the band that is most dispersive in the 
z direction. Bands 5 through 8 have predominantly C p z 
origin (band 5 with a considerable hybridization with Yb 
/), while bands 1 and 3 show only small contributions 
within either C or Yb muffin-tin spheres, but rather in 
the interstitial areas, thus being mostly free-electron like. 




the lower index is blue (darker in grey scale), (color online) 



While at the Fermi level most of the DOS (inset in 
Fig. 3) comes from the extended states in the inter- 
stitial region, there is a noticeable share of the Yb / 
states, which rapidly grows with increasing BE. Contri- 
bution of the Yb d states (mostly d z 2_ r 2 symmetry) is 
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also considerable. This indicates that (i) Yb phonons 
may play a role in superconductivity and (ii) substitut- 
ing Yb by other divalent elements, like Ca, may change 
superconducting properties substantially, or destroy su- 
perconductivity at all. The total density of states at the 
Fermi level is N(E F ) = 1.7 states/ (eV-formula). The av- 
erage Fermi velocity in- and out-of-plane, respectively, is 
4.8xl0 7 and 2.7xl0 7 cm/sec. The corresponding com- 
ponents of the plasma frequency are 6.2 and 3.6 eV, im- 
plying, in the constant scattering rate approximation, a 
resistivity anisotropy of about 3. 

The large unit cell of YbCg allows for multiple zone- 
center modes, which in principle can be computed by the 
frozen phonon method. 24 modes corresponding to in- 
plane motion of C are expected to be very hard, since 
C p x>y orbitals are removed from the Fermi level (as op- 
posed, for instance, to MgB 2 ). Indeed, we calculated the 
frequency of the fully symmetric A lg mode to be 1540 
cm -1 . Most of the 12 out-of plane C modes, by sym- 
metry, have zero deformation potential at T and thus 
do not couple with electrons, while the E g mode low- 
ers the symmetry to triclinic and is prohibitively diffi- 
cult to compute accurately. Thus, we are left with the 
6 Yb modes, of which the three even modes correspond 
to a double degenerate Ei g phonon involving displace- 
ments of Yb parallel to the graphite plane, and a B\ g 
one, corresponding to Yb displacements along c. These 
modes are rather soft; our calculations place them at 77 
and 153 cm -1 , respectively. Unfortunately, the electron- 
phonon coupling exactly at the zone center is strongly 
suppressed by symmetry for these three modes (if the al- 
ternating Yb planes were not shifted with respect to each 
other, cf. Fig. 1, it would be forbidden, so in the actual 
structure it is suppressed to the extent of the weakness 
of Yb-Yb interaction. However, their low frequency and 
considerable presence of Yb character at the Fermi level 
suggest that coupling of these and similar modes, when 
integrated over the Brillouin zone, may be sufficient to 
explain superconductivity in YbC6- 
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